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The central inclusive jet cross section has been measured using a successive-combination algorithm 
for reconstruction of jets. The measurement uses 87.3 pb _1 of data collected with the D0 detector 
at the Fermilab Tevatron pp Collider during 1994-1995. The cross section, reported as a function 
of transverse momentum (p T > 60 GeV) in the central region of pseudorapidity (\rj\ < 0.5), exhibits 
reasonable agreement with next-to-leading order QCD predictions, except at low pr where the 
agreement is marginal. 
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Jet production in hadronic collisions is understood 
within the framework of Quantum Chromodynamics 
(QCD) as a hard scattering of constituent partons 
(quarks and gluons), that, having undergone the interac- 
tion, manifest themselves as showers of collimated parti- 
cles called jets. Jet finding algorithms associate clusters 
of these particles into jets so that the kinematic prop- 
erties of the hard-scattered partons can be inferred and 
thereby compared to predictions from perturbative QCD 
(pQCD). 

Historically, only cone algorithms have been used to 
reconstruct jets at hadron colliders Although well- 
suited to the understanding of the experimental system- 
atics present in the complex environment of hadron col- 
liders, the cone algorithms used in previous measure- 
ments by the Fermilab Tevatron experiments ^,|| present 
several difficulties: an arbitrary procedure must be imple- 
mented to split and merge overlapping calorimeter cones, 
an ad- hoc parameter, lZ scp [H, is required to accommo- 
date the differences between jet definitions at the parton 
and detector levels, and improved theoretical predictions 
calculated at the next-to- next-to-leading-order (NNLO) 
in pQCD are not infrared safe, but exhibit sensitivity to 
soft radiation [||. 

A second class of jet algorithms, which does not suf- 
fer from these shortcomings, has been developed by sev- 
eral groups These recombination algorithms suc- 
cessively merge pairs of nearby objects (partons, parti- 
cles, or calorimeter towers) in order of increasing relative 
transverse momentum. A single parameter, D, which ap- 
proximately characterizes the size of the resulting jets, 
determines when this merging stops. No splitting or 
merging is involved because each object is uniquely as- 
signed to a jet. There is no need to introduce any ad-hoc 
parameters, because the same algorithm is applied at the 
theoretical and experimental level. Furthermore, by de- 
sign, clustering algorithms are infrared and collinear safe 
to all orders of calculation. In this Letter, we present the 
first measurement of the inclusive jet cross section using 
the k± algorithm |]-[|] to reconstruct jets at the y/s = 1.8 
TeV Tevatron proton-antiproton collider. 

The differential jet cross section was measured in bins 
of pt and pseudorapidity, r\ = — ln[tan(0/2)], where 
6 is the polar angle relative to the z axis pointing in 
the proton beam direction. The k± algorithm imple- 
mented at D0 Q is based on the clustering algorithm 
suggested in Ref. The algorithm starts with a list of 
objects. For each object with transverse momentum PT,h 
we define da = p T ii an d for each pair of objects, rfy = 
mm(p Ti ,p T j) (ARij) 2 /D 2 , where D is the free parame- 
ter of the algorithm and (ARij) 2 — (Acj)^) 2 + (Arnj) 2 is 
the square of their angular separation. If the minimum of 
all da and dij is a dij, then the objects i and j are com- 
bined, becoming the merged four- vector (Ei+Ej, pi+Pj). 
If the minimum is a da, the object i is defined as a jet 



and removed from subsequent iterations. This procedure 
is repeated until all objects are combined into jets. Thus 
k± jets do not have to include all objects in a cone of 
radius D, and they may include objects outside of this 
cone. 

The primary tool for jet detection at D0 is the liquid- 
argon/uranium calorimeter JTofl , which has nearly full 
solid-angle coverage for \r]\ < 4.1. The first stage (hard- 
ware) trigger selected inelastic collisions as defined by sig- 
nal coincidence in the hodoscopes located near the beam 
axis on both sides of the interaction region. The next 
stage required energy deposition in any Ai] x A</> = 0.8 x 
1.6 region of the calorimeter corresponding to a trans- 
verse energy (Er) above a preset threshold. Selected 
events were digitized and sent to an array of processors. 
At this stage jet candidates were reconstructed with a 
cone algorithm (with radius R = [(Ar;) 2 + (A0) 2 ]2 = 
0.7), and the event was recorded if any jet Er exceeded 
a specified threshold. Jet Er thresholds of 30, 50, 85, 
and 115 GeV accumulated integrated luminosities of 0.34, 
4.46, 51.5, and 87.3 pb _1 , respectively 0. 

Jets were reconstructed offline using the k± algorithm, 
with D — 1.0. This value of D was chosen because, 
at next-to-leading-order (NLO), it produces a theoreti- 
cal cross section that is essentially identical to the cone 
prediction for R = 0.7 ||], which D0 used in its previ- 
ous publications on jet production ||. The vertices of 
the events were reconstructed using the central tracking 
system JhJ . A significant portion of the data was taken 
at high instantaneous luminosity, where more than one 
interaction per beam crossing was probable. When an 
event had more than one reconstructed vertex, the quan- 
tity St = |Sp*p et | was defined for the two vertices that 
had the largest numbers of associated tracks, and the ver- 
tex with the smallest St was used for calculating all kine- 
matic variables pjlljl- To preserve the pseudo-pro jective 
nature of the D0 calorimeter, the vertex z-position was 
required to be within 50 cm of the center of the detector. 
This requirement rejected (10.6 ± 0.1)% of the events, 
independent of jet transverse momentum. 

Isolated noisy calorimeter cells were suppressed with 
online and offline algorithms fuJl . Background intro- 
duced by electrons, photons, detector noise, and accel- 
erator losses that mimicked jets were eliminated with jet 
quality cuts. The efficiency of the jet selection is approx- 
imately 99.5% and nearly independent of jet pt- The 
imbalance in transverse momentum, "missing transverse 
energy," was calculated from the vector sum of the E XjV 
values in all cells of the calorimeter. Background from 
cosmic rays or incorrectly vertexed events was eliminated 
by requiring the missing transverse energy in each event 
to be less than 70% of the pt of the leading jet. This 
criterion caused essentially no loss in efficiency. 

The D0 jet momentum calibration H, applied on a 
jet-by-jet basis, corrects on average the reconstructed 
calorimeter jet pTto that of the final-state particles in the 
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FIG. 1. The central (|?7| < 0.5) inclusive jet cross section 
obtained with the k± algorithm at y/s = 1.8 TeV. Only sta- 
tistical errors are included. The solid line shows a prediction 
from NLO pQCD. 

jet. The correction accounts for contribution from back- 
ground from spectator partons (the "underlying-event," 
determined from minimum-bias events), additional inter- 
actions, pileup from previous pp crossings, noise from 
uranium radioactivity, detector non-uniformities, and for 
the global response of the detector to hadronic jets. Un- 
like the cone algorithm, the k± algorithm does not re- 
quire additional corrections for showering in the calorime- 
ter ||. For \rj\ < 0.5, the mean total multiplicative 
correction factor to an observed pt of 100 (400) GeV is 
1.094 ±0.015 (1.067 ± 0.020). 

The inclusive jet cross section for |r;| < 0.5 was cal- 
culated over four ranges of transverse momentum, each 
using data from only a single trigger threshold. The more 
restrictive trigger was used as soon as it became fully ef- 
ficient. The average differential cross section for each 
Pt bin, ePer / '(dpTdrj), was measured as N/(Ai]ApTeL), 
where Ar) and Apx are the rj and pr bin sizes, N is the 
number of jets observed in that bin, e is the overall effi- 
ciency for jet and event selection, and L represents the 
integrated luminosity of the data sample. 

The measured cross section is distorted in px by the 
momentum resolution of the D0 calorimeter. The frac- 
tional momentum resolution was determined from the 
imbalance in pr in two-jet events jy]. Although the res- 
olution in jet pt is essentially Gaussian, the steepness 
of the pt spectrum shifts the observed cross section to 
larger values. At 100 (400) GeV, the fractional resolution 
is 0.061±0.006 (0.039±0.003). The distortion in the cross 
section due to the resolution was corrected by assuming 
an ansatz function, Ap^. B (l — 2pT/y/s) c , smearing it 
with the measured resolution, and fitting the parameters 
A, B and C so as to best describe the observed cross 
section. The bin-to-bin ratio of the original ansatz to the 
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TABLE I. Inclusive jet cross section of jets reconstructed 
using the k± algorithm in the central pseudorapidity region 
(M<0.5). 



smeared one was used to remove the distortion due to 
resolution. The unsmearing correction reduces the ob- 
served cross section by (5.7 ± 1)% ((6.1 ± 1)%) at 100 
(400) GeV. 

The final, fully corrected cross section for |»;| < 0.5 is 
shown in Fig. [l], along with the statistical uncertainties. 
Listed in Table || are the pr range, the best pr bin cen- 
troid, the cross section, and uncertainties in each bin. 
The systematic uncertainties include contributions from 
jet and event selection, unsmearing, luminosity, and the 
uncertainty in the momentum scale, which dominates at 
all transverse momenta. The fractional uncertainties for 
the different components are plotted in Fig. || as a func- 
tion of the jet transverse momentum. 

The results are compared to the pQCD NLO prediction 
from jetrad | fi2[ , with the renormalization and factor- 
ization scales set to p™ ax /2, where p™ ax refers to the proi 
the leading jet in an event. The comparisons are made 
using parametrizations of the parton distribution func- 
tions (PDFs) of the CTEQ ftf] and MRST (l| families. 
Figure || shows the ratios of (data-theory) /theory. The 
predictions lie below the data by about 50% at the lowest 
p T and by (10 - 20)% for p T > 200 GeV. To quantify the 
comparison in Fig. |3|, the fractional systematic uncertain- 
ties are multiplied by the predicted cross section, and a 
X 2 comparison, using the full correlation matrix, is car- 
ried out Q . The results are shown in Table ||. Though 
the agreement is reasonable (% 2 /dof ranges from 1.56 to 
1.12, the probabilities from 4 to 31%), the differences in 
normalization and shape, especially at low pr, are quite 
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FIG. 2. Fractional experimental uncertainties on the cross 
section. The discontinuities in the luminosity uncertainty are 
related to the use of different triggers [11]. 



PDF 


x 2 


X 2 /dof 


Probability (%) 


MRST 


26.8 


1.12 


31 


MRSTgT 


33.1 


1.38 
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MRSTgX 


28.2 


1.17 
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CTEQ3M 


37.5 


1.56 
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CTEQ4M 


31.2 


1.30 


15 


CTEQ4H.T 


27.2 


1.13 


29 


MRST+Hadroniz. 


24.0 


1.00 


46 


CTEQ4HJ+Hadroniz. 


24.3 


1.01 


44 



2 comparison (24 degrees of freedom) between 



TABLE II. x 

jetrad, with renormalization and factorization scales set to 
p™ ax /2, and data for various PDFs. The last entries include 
a hadronization correction obtained from herwig (see text). 



large. The points at low pt have the highest impact on 
the x 2 . If the first four data points are not used in the \ 2 
comparison, the probability increases from 29% to 77% 
when using the CTEQ4HJ PDF. 

While the NLO predictions for the inclusive cross sec- 
tion for k± (D = 1.0) and cone jets (R — 0.7, lZ scp = 1.3 
in the same \r]\ < 0.5 interval are within 1% of each other 
for the pt range of this analysis JTTJ , the measured cross 
section using k± is 37% (16%) higher than the previously 
reported cross section using the cone algorithm jl]| at 60 
(200) GeV. This difference in the cross sections is con- 
sistent with the measured difference in pr for cone jets 
matched in rj — <j> space to k± jets. k± jets were found to 
encompass 7% (3%) more transverse energy at 60 (200) 
GeV than cone jets |Q. 

The effect of final-state hadronization on reconstructed 
energy, which might account for the discrepancy between 
the observed cross section using k± and the NLO pre- 
dictions at low pt, and also for the difference between 
the k± and cone results, was studied using herwig (ver- 
sion 5.9) (l6) simulations. Figure [| shows the ratio of 
Pt spectra for particle-level to parton-level jets, for both 
the kj_ and cone algorithms. Particle cone jets, recon- 
structed from final state particles (after hadronization), 
have less prthan the parton jets (before hadronization), 
because of energy loss outside the cone. In contrast, 
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FIG. 3. Difference between data and jetrad pQCD, nor- 
malized to the predictions. The shaded bands represent the 
total systematic uncertainty. In the bottom plot a herwig 
hadronization contribution has been added to the prediction 
(open circles). 



particle jets are more energetic than their progenitors 
at the parton level, due to the merging of nearby partons 
into a single particle jet. Including the hadronization 
effect derived from herwig in the NLO jetrad predic- 
tion improves the x 2 probability from 29% to 44% (31% 
to 46%) when using the CTEQ4HJ (MRST) PDF. Wc 
have also investigated the sensitivity of the measurement 
to the modeling of the background from spectator par- 
tons through the use of minimum bias events, and found 
that it has a small effect on the cross section: at low pt, 
where the sensitivity is the largest, an increase of as much 
as 50% in the underlying event correction decreases the 
cross section by less than 6%. 

In conclusion, we have presented the first measurement 
in proton- antiproton collisions at y/s — 1.8 TeV of the 
inclusive jet cross section using the k± algorithm. Quan- 
titative tests show reasonable agreement between data 
and NLO pQCD predictions, except at low pt where the 
agreement is marginal. The degree of agreement can be 
slightly improved by incorporating a hadronization con- 
tribution of the kind predicted by herwig. 
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